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ABSTRACT 



Context. Young stars are born with envelopes, which in the early stages obscure the central (proto)star and circumstellar disk. In the 
Class I stage, the disks are still young, but the envelopes are largely dispersed. This makes the Class I sources ideal targets for studies 
of the early stages of disks. 

Aims. We aim to determine the masses of the envelopes, disks, and central stars of young stellar objects (YSOs) in the Class I stage. 
Methods. We observed the embedded Class I objects IRS 63 and Elias 29 in the p Ophiuchi star-forming region with the 
Submillimeter Array (SMA) at 1 . 1 mm. 

Results. IRS 63 and Elias 29 are both clearly detected in the continuum, with peak fluxes of 459 and 47 mjy/beam, respectively. 
The continuum emission toward Elias 29 is clearly resolved, whereas IRS 63 is consistent with a point source down to a scale of 3" 
(400 AU). The SMA data are combined with single-dish data, and both disk masses of 0.055 and < 0.007 M Q and envelope masses 
of 0.058 and < 0.058 M Q are empirically determined for IRS 63 and Elias 29, respectively. The disk+envelope systems are modelled 
with the axisymmetric radiative-transfer code RADMC, yielding disk and envelope masses that differ from the empirical results by 
factors of a few. HCO + J = 3-2 is detected toward both sources, HCN J = 3-2 is not. The HCO + position-velocity diagrams are 
indicative of Keplerian rotation and allow an estimate of the mass of the central stars. For a fiducial inclination of 30°, we find stellar 
masses of 0.37 ± 0.13 for IRS 63 and 2.5 ± 0.6 M G for Elias 29. 

Conclusions. The sensitivity and spatial resolution of the SMA at 1.1 mm allow a good separation of the disks around Class I YSOs 
from their circumstellar envelopes and environments, and the spectral resolution makes it possible to resolve their dynamical structure 
and estimate the masses of the central stars. The ratios of the envelope and disk masses M cnv /M disk are found to be 0.2 for IRS 63 and 
6 for Elias 29. This is lower than the values for Class sources, which have M cnv /A/disk > 10, suggesting that this ratio is a tracer of 
the evolutionary stage of a YSO. 

Key words, circumstellar matter - planetary systems: protoplanetary disks - stars: formation - stars: individual (IRS 63, Elias 29) 



1. Introduction 

Low- and intermediate-mass stars are formed from the gravia- 
tional collapse of molecular cloud cores. In the earliest stages, 
the newly-formed protostar remains embedded in the remnants 
of this core, a cold envelope of dust and gas, which is gradually 
accreted by the young star (e.g., Shu 1977). Due to the angu- 
lar momentum initially present in the core, most of the enve- 
lope material does not fall directly onto the central protostar but 
is piled-up in a circumstellar disk (e.g., Terebey et al. 1984). 
Understanding this interplay between star, disk, and envelope is 
crucial in order to be able to relate the initial conditions of star 
formation such as the mass of the protostellar core to the end- 
product - namely the properties of the young star, and the mass 
and thus potential of the disk for forming planets. Some of the 
key questions include: Where does most of the mass reside at a 
given time? Will all the mass that is seen in prestellar cores or 
in the envelopes around deeply embedded sources end up in the 
star, or will a large fraction be dispersed from the system? How 
do the masses of the circumstellar envelope, the disk, and that of 
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the central star evolve over time, and how long does it take for 
the circumstellar matter to be accreted onto the star? 

Young stellar objects (YSOs) are usually classified according 
to their slopes in the infrared (IR) wavelength regime. Originally, 
the LW classification ( |Lada & Wilking||1984"| |Lada||1987| > ran 
from the embedded Class I, via the optically visible Class II or 
classical T Tauri stars, to the Class III spectral energy distribu- 
tions of post-T Tauri stars. Later, the Class stage was added to 
this classification (see, e.g., |Andre et al. 1993| l, where the Class 
sources are distinguished from the Class I sources through their 
high relative luminosity at submillimetre (submm) wavelengths. 
This classification roughly reflects the evolutionary stage of the 
YSOs under consideration. The most deeply embedded Class 
sources are thought to evolve through the Class I stage while dis- 
sipating their circumstellar envelopes. Eventually they become 
optically visible as pre-main sequence T Tauri stars with circum- 
stellar disks. 

The Class and Class II YSOs have been studied quite ex- 
tensively with high-resolution (sub)millimetre interferometers 
(e.g., |j0rgensen et al.|2007||Andrews & Wil liams 2007 and ref- 
erences therein). Studies of deeply embedded Class YSOs have 
shown that circumstellar disks are formed early (Harvey et al. 
|2003 ; J0rgens en et al.|200 4), but in these systems, it is difficult 
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to separate the emission from the disk from that of the envelope. 
The Class I sources, in which a large part of the original circum- 
stellar envelopes has been dissipated, are ideal to study young 
disks in YSOs. Interferometric studies of Class I objects have so 
far been largely limited to studies at around 3 mm (e.g.,[Ohashi 



etal.|1997| [Hogerheijde 'et al.|1997| [T9981 |Looney et al.]2 000) 
The Submillimeter Array (SMA) allows observations around 
1 mm, where the thermal dust continuum emission is an order 
of magnitude stronger than at 3 mm. Also at these shorter wave- 
lengths it is possible to detect the higher rotational transitions 
of the molecules that trace the dense gas in the disks and inner 
envelopes of the young systems, rather than lower-density ex- 
tended envelope emission. 

We here present SMA observations of two Class I objects 
that appear to be in an evolutionary stage where the kinematics 
of the circumstellar material are no longer dominated by infall. 
IRS 63 (WLY 2-63, GWAYL 4) and Elias 29 (Elia 2-29, WLY 
1-7) are located in the p Ophiuchi cloud, taken to be at a dis- 
tance D = 125 + 25 pc ( |de Geus et al. 1989 1. Spitzer photometry 
from the "Cores to Disks" Legacy program ( |Evans et al.|2 003 ) 
was used to determine the source's infrared colours and their 
bolometric luminosities and temperatures. IRS 63 and Elias 29 
have similar bolometric temperatures, 7boi = 351 K and 391 K, 
respectively, which places them in the LW Class I regime. The 
bolometric luminosities were calculated to be Lboi = 0.79 L Q 
for IRS 63 and 13.6 L Q for Elias 29. According to their infrared 
slopes of a2-24iim =0.15 and 0.42, IRS 63 is the slightly more 
evolved, falling into the "flat-spectrum" class that separates the 
Class I from the Class II sources (Gr eene et al.|1994] ). The val- 
ues for a iR quoted here are from un-dereddened observations. 
Correction for the extinction toward the p Ophiuchi star-forming 
region (e.g., |Flaherty et al.|2007| l would result in lower values for 
am; in other words, the sources may be slightly more evolved 
than the raw values of o^r suggest. In terms of environment, 
the two sources are quite each other's opposites. The SCUBA 
850 /jm map (| Johnstone et al.|2004| > from the COMPLETE sur- 
vey ( |Ridge et a l. 2006 ) shows that IRS 63 is an isolated compact 
source. Elias 29, on the other hand, is located in a dense ridge 
of molecular material, which contains several more YSOs. It is 
likely that these YSOs were formed from condensations in this 
molecular ridge. 

We present the masses of all main components of Class I 
YSOs, i.e., the central star, the disk, and the envelope, for the first 
time. IRS 63 and Elias 29 are the first two sources in a larger sur- 
vey of Class I objects studied with the SMA at 1.1 mm. These 
observations complement the survey of Class sources in the 
PROS AC programme (J0rgensen et al. 2007) and will allow us 
to trace the similarities and differences of these evolving proto- 
stars. The results of the complete campaign will be presented in 
a future paper. In £j2] the observations are presented, and the re- 
sults and implications are discussed in §f|3]and|4] We summarise 
the main conclusions in Sj5] 

2. Observations 



IRS 6 3 and Elias 29 were observed with the SMAQ ( |Ho et al. 
2004| l on 15 and 17 May 2006, respectively. Weather conditions 
were good on 15 May, with zenith optical depths at 225 GHz 
T225 = 0.04-0.06 (as measured by a tipping radiometer operated 



by the Caltech Submillimeter Observatory). Conditions were 
slightly worse on 17 May, with T225 starting at 0.16, falling to 0.1 
during the first half of the night. Physical baselines ranged from 
1 1.6 to 62.0 metres, and the resulting projected baselines ranged 
from 5.9 to 63.8 LI, yielding a resolution of about 4.0 x 2.3" 
(natural weighting). The correlator was configured to observe 
the lines of HCO + J = 3-2 and HCN J = 3-2 together with the 
continuum emission at 268 and 278 GHz. The line data were 
taken with 512 channels of 0.2 MHz width each, resulting in a 
velocity resolution of 0.23 km s . 

Calibration was done using the MIR package (Qi 2005 [l, 
and imaging was done using the MIRIAD package ( |Sault et Ik] 
[1995]). The quasars QSO B1622-297 and QSO B1514-241 
served as gain calibrators, and the absolute fluxes were cali- 
brated on Uranus. The absolute flux calibration is estimated to 
have an uncertainty of ~20 %. The passbands were calibrated on 
Callisto and the quasar QSO B1253-055 (3C 279). 

3. Results 

The basic results of the observations are shown in Fig. [T] and 
summarised in Table Q] 



3.1. Continuum data 

Both sources are clearly detected in the continuum at 1.1 mm. 
Point-source and Gaussian fits were done in the (m, v) plane to 
determine the peak and integrated fluxes, respectively. A con- 
tinuum peak flux of 459 mJy bm _1 and an integrated flux of 
474 mJy were found for IRS 63, at a position which agrees with 
the 2MASS K s -band po sition within 071. IRS 63 wa s also ob- 
served with the SMA by | Andrews & Williams] ( |2007| ) at higher 
resolution at 1.3 mm. Their interferometer flux for IRS 63, to 
which they refer as LI 709 B^] is consistent with the value here 
and a mm slope a mm = 2.0 + 1.0. The data of Elias 29 are best 
fitted by two sources. The brightest of the two agrees with the 
2MASS Ks-band position of Elias 29 within 072, whereas the 
other source is offset by 3" . The second peak is attributed to an 
enhancement in the ridge, from which Elias 29 likely formed, 
and is designated by "Ridge" in Table [T] Elias 29 and the ridge 
component yield continuum peak fluxes of 47 and 31 mJy bm , 
and integrated fluxes of 72 and 28 mJy, respectively. 

Plots of the visibility amplitudes as functions of projected 
baseline for IRS 63 and Elias 29 are shown in Fig. [2] Elias 29 
shows a steeply rising flux toward the shortest baselines, indicat- 
ing that an envelope is present around this source. The emission 
toward IRS 63 does not change appreciably with baseline length, 
indicating that this source is unresolved up to baselines of 60 kA, 
or down to physical scales of about 400 AU. 

Fig. [2] also shows the estimated single-dish flux at 1.1 mm. 
For this, the "mapped" fluxes from Andre & Montmerle (1994) 
at 1.25 mm were combined with newly determined 850 jum 
fluxes (SCUBA maps from the COMPLETE survey, |Johnstone| 
et al.|2004| |Ridge et al.|200"6] ). The SEDs are assumed to have 
a F v oc V dependence in this wavelength range, and are in- 
terpolated to the 1.1 mm wavelength at which the SMA ob- 
servations were conducted. Fluxes of 355 and 780 mJy are 
found for Elias 29 and IRS 63, respectively, with an estimated 



1 The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 



2 IRS 63 and L1709 B are often referred to as the same source (e.g., 
Andre & Montmerle 1994). However, L1709 B would appear to lie 
about 4 arcmin north of IRS 63 (e.g. [Ben son & Myers 1989). The 
source that I Andrews & Williams] (20071 ) refer to as LI 709 B, is in fact 
IRS 63. 
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Fig. 1. SMA images of the 1.1 mm continuum (top panels) and 
integrated HCO + 3-2 (bottom panels) emission. Contours are 
drawn at 2, 4, 8, 16, 32, and 64 times the respective RMS noise 
levels, with the RMS being about 6 mJy birr 1 for the continuum 
and about 65 mJy birr 1 for the line data. Negative contours are 
dashed. The positional offsets are with respect to the coordinates 
of the observational phase centre. The dashed lines indicate how 
the position velocity diagrams (Fig. |7]i are directed. 
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(u s +v 2 ) 1/2 [kA] 

Fig. 2. Continuum amplitude as a function of projected baseline 
for IRS 63 (left panel) and Elias 29 (right panel). The data points 
give the amplitude per bin, where the data are binned in annuli 
according to (u, v) distance. The error bars show the statistical 
lcr errors, most often smaller than the data points, and the dot- 
ted lines give the expected amplitude for zero signal, i.e., the 
anticipated amplitude in the absence of source emission. The 
half-open circles at zero (u, v) distance give the zero-spacing 
1 . 1 mm flux , interpolated between 850 pm and 1 .25 mm single- 
dish fluxes ([Andre" & Montmerle||1994[ |Johnstone et al.| [2Q04 ; 
Ridge et al.|2006[ see text). 
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Fig. 3. First-moment map of IRS 63. Contours are drawn at 2, 4, 
and 8 times the RMS of 67 mJy birr 1 . The positional offsets are 
with respect to the coordinates of the observational phase centre, 
and the cross indicates the continuum position of IRS 63 from a 
point-source fit in the (u, v) plane, see Table [T] A colour version 
of this figure can be obtained from the electronic version of the 
paper. 

3.2. Line data 

HCO + J = 3-2 is detected toward both sources, whereas HCN J 
= 3-2 is not. Elias 29 appears to be about ten times as strong in 
HCO + J = 3-2 integrated emission as IRS 63 (see Table[T]i. The 
emission toward IRS 63 is compact, coincident with the con- 
tinuum peak, and shows a velocity gradient in the north-south 
direction, see Fig. [3] The integrated emission in the direction 
of Elias 29 is more extended and can be fitted with two el- 
liptical Gaussians, as is shown in the left panel of Fig. [4] The 
centre position of the smaller Gaussian coincides with the con- 
tinuum peak, as well as with the infrared position, and is at- 
tributed to Elias 29 itself. The second Gaussian has a very elon- 
gated shape, offset from the infrared source, and is attributed 
to a density enhancement in the nearby ridge. The right panel 
of Fig. [4] shows the HCO + J = 3-2 spectra toward the contin- 
uum positions of Elias 29 and IRS 63, binned to 0.9 km s _1 . The 
HCO + J = 3-2 emission toward IRS 63 appears to be centred at 
Vlsr = 3.3 kms -1 . 

4. Discussion and interpretations 

4. 1. Envelope and disk masses 

The basic data as presented in the previous section are used 
to determine the masses of the disks and the envelopes in the 
both empirically (Sect 



systems, both empirically (Sect. 4.1.1 
radiative-transer modelling (Sect. 4.1.2 
tween the two methods and the over-a 
cussed in Sect. 14.1.31 



and through detailed 
1. The differences be- 
1 conclusions are dis- 



uncertainty of 25%. These fluxes are significantly larger than 
the interferometer fluxes and indicate that extended emission 
is present around the compact components picked up with the 
SMA. This extended emission may be due to a circumstellar en- 
velope, or to surrounding or intervening interstellar clouds, since 
the "mapped" regions of Andre & Montmerle (1994) are rather 
large, 40" for Elias 29 and 60" for IRS 63. 



4.1.1. Empirical results 

A first-order estimate of the disk and envelope masses can be 
obtained by assuming that the continuum flux of the envelope 
is resolved out by the interferometer on the longest baselines 
(Fig. |2j. Hence, the emission on the longest baselines is solely 
due to the disk, whereas the single-dish flux comes from the disk 
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Fig. 4. Left panel: HCO + J = 3-2 integrated emission toward 
Elias 29. The grayscale is linear from 0.2 to 1.2 Jy birr 1 . The 
emission is fitted with two Gaussians, the one in the north-south 
direction tracing the emission due to Elias 29, and the other, 
larger one tracing emission from a density enhancement in the 
ridge. Right panel: spectral lines, integrated over 2 x 2" regions, 
centred on the continuum positions of Elias 29, of the ridge 
(shifted by +4 Jy bor 1 ), and of IRS 63 (shifted by -4 Jy bnr 1 ), 
and binned to 0.9 km s . The HCO + J = 3-2 emission toward 
IRS 63 appears to be centred at Vlsr = 3.3 km s . 

Table 1. Results of SMA observations at 1.1 mm. 





IRS 63 


Elias 29 


Ridge 


F„ (P) a [mJy birr 1 ] 


459 


47 


31 


RA a [J2000] 


16:31:35.65 


16:27:09.42 


16:27:09.48 


Dec a [J2000] 


-24:01:29.56 


-24:37:18.91 


-24:37:22.48 


RMS b [mJy bm-'] 


6 


5 


5 


Circular Gaussian fits to all data 


F v (Gf [mJy] 


474 


72 


28 


FWHM C [arcsec] 


0.55 ± 0.08 


2.2 ±0.5 




Circular Gaussian fits to data >16 kA (scales < 13 arcsec) 


F v (Gf [mJy] 


480 


57 


27 


FWHM C [arcsec] 


0.62 ± 0.08 


1.7 ±0.7 




Line data 


HCO +d [Kkms 1 ] 


3.1 ±0.3 


23.7 ±2.5 


23.5 ±2.5 


HCN d [Kkms- 1 ] 


< 1.8 


< 3.8 


< 3.8 



Point source fit in the (u, v) plane. 
Calculated from the cleaned image. 
Circular Gaussian fit in the («, v) plane. 

Integrated intensities are from a 4 x 4" square around the contin- 
uum emission. The synthesised beam is ~ 4 x 2.5". HCN was not 
detected; quoted values are 3cr upper limits. 



and envelope combined. Assuming that the disk and the envelope 
are isothermal and the dust emission is optically thin, the disk 
mass is given by 



Ms 



Fy'VD 2 

K v B v (7dust) 



(1) 



where F v is the flux at frequency v on the longest baselines, 
is the gas-to-dust ratio, k v is the dust opacity at v, and B v (T& nsl ) 
is the emission from a black body at Ti ust . Using fiducial val- 
ues for all these parameters - T* = 100, Tjust = 30 K, and 
k v = 1.18 cm 2 g -1 ("OH5" coagulated dust with icy mantles, 
found in the fifth column of Table 1 of Ossenkopf & Henning 
1994| interpolated to the observed frequency v = 273 GHz) - 
disk masses of 0.055 M G and 0.007 Mq are found for IRS 63 
and Elias 29, respectiv ely. Using the opacity law as quoted by 

(1990i, k v = 10{v/l0 12 Hzf, and taking £ = 1, 



a dust opacity at 273 GHz of k v = 2.73 cm 2 g 1 is found, and 
the disk masses decrease to 0.024 M G for IRS 63 and 0.003 M G 
for Elias 290Note that this method of estimating the disk mass 
is only valid when the disk is unresolved by the observations. 
If the disk were clearly resolved, the flux on the longest base- 
lines would not include the emission from the outer disk, but 
only that from smaller radii, unresolved by those baselines. Due 
to its closeness, the enhancement in the ridge near Elias 29 may 
contribute to its flux, even on the longest baselines of > 50 k/l or 
scales < 4". However, in the estimate of the flux on the longest 
baselines (Table [TJ, the ridge component was fitted simultane- 
ously with Elias 29, and hence its contribution to the derived 
disk mass is expected to be minimal. 

In a similar fashion the envelope masses in the systems 
can be found. First, the contributions from the envelopes to the 
1.1 mm single-dish fluxes are determined by subtracting the 
fluxes at baselines > 16 k/l (see Table [T} from the single-dish 
fluxes as found in Sect. [3] The envelope masses are then found 
by using Eq. (TfJ), with the temperatures of the envelopes taken to 
be 20 K. For an optically thin envelope around stars with lumi- 
nosities of 0.79 and 13.6 L Q applicable to IRS 63 and Elias 29, 
a temperature of 17-29 K is expected at a radius of 940 AU 
corresponding to the size of the JCMT single-dish beam (from, 
e.g., Eq. (2) of Chandler & Richer 2000 taking the opacity in- 
dex /3 — 1). This yields an envelope mass M env = 0.058 Mo for 
IRS 63 within a 30" radius, where the contribution from the en- 
velope to the total 1.1 mm flux is 38%. Coincidentally, also for 
Elias 29 an envelope mass M env = 0.058 M G , within a 20" ra- 
dius, is found. However, in this system the contribution from the 
envelope to the total 1.1 mm flux is 84%. As for the disk mass, 
also the envelope mass found for Elias 29 should be taken as an 
upper limit because of the closeness of the ridge enhancement. 
Note, however, that it is unlikely that the bulk of the continuum 
emission is due to this enhancement, because of the coincidence 
of the continuum peak with the infrared source and also with 
the single-dish submillimetre positions. Using the opacities of 
Beckwith et al. ( 1990) ) will lower all masses by a factor of about 
2. 

4.1.2. Modelling 

To provide more physical grounding to the values for Mdi S k 
and M env estimated in the previous section, a radiative transfer 
code was used to fit simultaneously the millimetre interferom- 
etry and the spectral energy distributions. IRS 63 and Elias 29 
were modelled with a new version of the programme RADMC 
(e.g., |Dullemond & Dominik|2004[ ). RADMC is an axisymmet- 
ric Monte-Carlo code for dust continuum radiative transfer in 
circumstellar disks and envelopes in which the stellar photons 
are traced in three dimensions. The code is based on the method 
of |Bjorkman & Wood| ( p00T] i. 

Th e density structure of the disk is given by (see Crapsi etaL] 
2007) 



Pdisk(r, &) 



y/(2n)H(r) 



exp 



' 1 


r cos(6) 


2\ 


2 

V 


[ H(r) _ 


) 



(2) 



where r is the radial distance from the central star and 6 is the 
angle from the axis of symmetry. The inner radius was fixed to 
0.1 AU for IRS 63, and to 0.25 AU for Elias 29, to account for 



The opacity law quoted by 



Beckwith et al. 



19901 is k v 



Beckwith et al. 



0.l(v/10 12 Hz) for the gas and dust combined, implicitly assuming a 
gas-to-dust ratio = 100, consistent with our work. 
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Table 2. Model parameters for IRS 63 and Elias 29. 





IRS 63 


Elias 29 


Modelling parameters and results 


Stellar luminosity (fixed) 


0.79 L 


13.6 L 


Disk inner radius (fixed) 


0.1 AU 


0.25 AU 


Disk outer radius, Rq 


100 AU 


200 AU 


Disk height at outer radius, H /R (fixed) 


0.17 


0.17 


Envelope outer radius (fixed) 


10,000 AU 


10,000 AU 


Centrifugal radius, R rot 


100 AU 


300 AU 


Inclination 


30° 


30° 


Disk mass, M disk 


0.13 M s 


0.004 M 


Envelope mass, M" nv 


0.022 Mo 


0.025 M 


Empirical results 


M c a nv [M ] 


0.058 


< 0.058 


M disk [Mo] 


0.055 


< 0.007 


AC [Mo] 


0.37 ±0.13 


2.5 + 0.6 



a The envelope mass is defined as the mass, not contained in the 
disk, within a single-dish radius, i.e., 30'73750 AU for IRS 63, and 
20"/2500 AU for Elias 29. 

b Assuming an inclination i = 30° (§ 4.2 for discussion). 



the higher luminosity of the star. At these radii, the temperature 
is of the order of the dust sublimation temperature. Note that the 
exact sublimation temperature, and hence the disk inner radius, 
depends on the exact dust species and local density. The vertical 
scale height of the disk is given by H(r) = r ■ Hq/Rq ■ (r/Ro) OA1 , 
wher e Rp was left free to vary, and Hq/Rq was arbitrarily fixed to 
0.17. |Chiang & Goldreichl ( |l997| l find H(r)/r oc r 2/1 . However, 
that estimate is based on grey opacities and constant surface 
density, and somewhat smaller flaring is usually more realistic. 
Hence, we chose H{r) oc r ■ r 017 . So was scaled to the disk mass, 
which was also left as a free parameter. The envelope density 



follows the equation for a rotating, infalling model (Ulrich 1976 
Crapsi et al.(2007] > 



Penvelope^, 8) — P0 



1.5 



1 + 



cos 6 
cos 0Q 



0.5 



Jv rot / 

COS 6 R mt 2 
H COS Gq 

2 cos Bo r 



(3) 



where 6q is the solution of the parabolic motion of an infalling 
particle, R rot is the centrifugal radius of the envelope, and po is 
the density in the equatorial plane at R mt . po was scaled to ac- 
comodate the total envelope mass, and R Iol , which can have a 
significant influence on the amplitude as a function of baseline 
length, was left free to vary. The outer radius of the envelope 
was fixed to 10,000 AU, where the temperature is similar to that 
of the ambient interstellar cloud. The parameters, as summarised 
in Table [2] give the best-fit models by eye; a full x 1 minimiza- 
tion is not warranted by the relatively low resolution and S/N 
of the data. Our models do not have outflow cavities in the en- 
velopes. Note that the actual physical structure of the envelope 
(centrifugal radius, outflow cavities, etc.) may considerably af- 
fect the mid-IR SED and thus the derived inclination from the 
models (e.g., Whit ney et aLp 003 ; J0rg ensen et al.|2005) . 

Two different dust species are used in these models, "hot" 
dust where the temperature is higher than 90 K, and "cold" 
dust where the temperature is below 90 K. The opacities used 
are "OH2" (coagulated dust without ice mantles, Ossenkopf & 
Henning 1994) for the hot dust, and "OH5" (coagulated dust 
with thin ice mantles) for the cold dust. The central stars were 
represented by sources with luminosities of 0.79 L for IRS 63 




Elias 29 




20 40 60 20 40 60 
(u 2 +v 2 ) 1/a [kAl 

Fig. 5. Continuum amplitude as a function of projected baseline 
for IRS 63 (left panel) and Elias 29 (right panel), with models 
overplotted. The dashed lines show models, consisting of an en- 
velope only. They fit the single-dish (sub)mm data well, but fall 
off too fast to longer projected baselines. The dash-dotted lines 
are from disk-only models. They fit the interferometric data well, 
but do not have enough flux on the shortest baselines. The solid 
lines are for the models that provide a good fit to the SED and to 
the spatial information provided by the interferometric observa- 
tions. The dotted lines indicate the expected amplitude for zero 
signal. 



and 13.6 L Q for Elias 29, which were taken to be equal to 
the bolometric luminosities of the sources, and the temperature 
structures in the disks and envelopes were then calculated by 
propagating the stellar photons through the systems. 

The results of the modelling are shown in Fig. [5] A model 
with only a 0.055 M Q disk reproduces the amplitude as a func- 
tion of baseline quite well for IRS 63, but fails to produce enough 
single-dish flux. A 0.35 Mq envelope-only model, on the other 
hand, produces enough single-dish flux in the (sub)mm regime, 
but falls off too quickly toward longer baselines. To allow for 
a comparison with the empirical model (Sect. |4. \.\\ the enve- 
lope mass was defined as the mass, present within the single- 
dish aperture (i.e., a 30" radius for IRS 63 and a 20" radius for 
Elias 29), and not contained in the disk. A model with both a 
disk and an envelope is necessary to explain the observations of 
IRS 63, and the best model consists of a disk of 0. 13 M Q and an 
envelope of 0.022 M . This model also fits the observed IRS 63 
SED well, if a foreground extinction of A v = 7 mag is assumed, 
see Fig. [6] The region in which IRS 63 is located has a very high 
Ay of about 24 mag as found in extinction maps produced by 



the Spitzer "Cores to Disks" Legacy program ( |Evans et al. 2003 
2007 ). These maps measure the large scale (~ 5') cloud extinc- 



tion on the basis of background stars, and therefore includes an 
extra contribution besides envelope extinction. 

The amplitude as a function of baseline for Elias 29 is best 
explained by a 0.004 M Q disk and a 0.025 Mq envelope. For this 
source, no further attempt to model the SED was made, because 
the close presence of the ridge enhancement makes a model that 
is axisymmetric on scales larger than several 100 AU unapplica- 
ble. 



4.1.3. Discussion 

Even though the disk and envelope masses derived from the em- 
pirical and detailed model results differ by up to a factor of a 
few, we are confident that they are accurate within this range. For 
example, a significantly larger disk in the IRS 63 model would 
show a very steep fall-off in the amplitude as a function of base- 



4 Delivery of data from the c2d Legacy Project: IRAC and MIPS 
(Pasadena, SSC,|Evans et al.|2007b. 
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Fig. 6. Spectral energy distribution of IRS 63 using the model 
parameters of Table 2 extincted by Ay=7 mag (solid line). The 
dashed line shows the model star, which has a luminosity of 
0.79 L , and is represented by a black-body with a temperature 



of 3 000 K. Overplotted are JHKL photometry ( |Haisch et aL 
|2002| , MIPS 24 and 70 fim fluxes (c2cQ>, JCMT 850 Jm 60" 
integrated flux (COMPLETE), SMA 1.10 mm peak flu x (this 
work; square), and IRAM 1 .25 mm 60" integrated flux ( | Andre 
|& Montmerle|1994| >. 



line length, which is not observed (Fig. [5]). Likewise, an enve- 
lope significantly more massive than 0.022 M Q found for IRS 63 
would obscure the central star so much that the near- and mid- 
IR flux would be severly underestimated (Fig. [6]). In the further 
discussion, we will use the disk and envelope masses as found 
by the detailed modelling, but note that the overall conclusions 
do not change if the empirical results are used. 

The envelope masses of IRS 63 and Elias 29 are in the same 
range as those presented by, e.g., Hogerheijde et al. (1997) for a 



sample of embedded Class I sources in the Taurus-Auriga star- 
forming region. However, the envelope masses of the deeply 
embedded Class objects are found to be considerably higher, 



> 1M ( 


J0rgensen et al. 2002 


Shirley et al.|2002 


Young et al. 


2003; Hatchell et al.|2007 1. While some fraction of sources may 



simply originate from lower mass cores, it is clear that all Class 
objects must pass through a stage with M env ~ 0.1 M on their 
way to the pre-main sequence stage. A larger sample is needed to 
address the question whether the mass that was initially present 
in the envelope will first accumulate in the disk, or whether it 
will pass through the disk and onto the star directly, leaving the 
disk in a steady state through the entire Class I stage. 

An interesting quantity for these YSOs is the ratio 
M env /Mdi s k, as this may be a direct measure for their evolutionary 
stage. The values for the two objects under consideration here - 
0.2 for IRS 63 and 6 for Elias 29 - are quite different, although 
the value for Elias 29 is rather uncertain due to the contribu- 
tion of the ridge enhancement to the continuum emission. From 
these simple arguments, Elias 29 is the less evolved, rather like 
the deeply embedded Class sources in the PROSAC sample 
( |j0rgensen et al.||2007[ >, which show M env /M disk > 10. IRS 63, 



on the other hand, is well on its way toward the optically visible 
Class II or T Tauri stage, which have M env ~ 0. This confirms the 
notion that the Lada classification is primarily a phenomenolog- 
ical one, not necessarily representing the actual physical stage 
the YSO is in (see, e.g., |Crapsi et al.|2007] >. 



4.2. Keplerian rotation and stellar masses 

Both IRS 63 and Elias 29 show a velocity gradient in the HCO + 
emission, which is interpreted as the rotation of a circumstel- 
lar disk. The mass of the central object can be estimated from 
position-velocity diagrams along the major axis of the disk, as 
indicated in Fig. [T] An elliptical Gaussian was fitted to each 
channel in the HCO + data in the (u, v) plane, using the MIRIAD 
task uvfit. This yielded a best-fit position per channel with cor- 
responding uncertainties in right ascension and declination, and 
consequently a declination offset from the phase centre in the 
physical plane. Note that the fitting was carried out in the (w, v) 
plane, so as not to be hindered by artefacts that may arise in the 
deconvolution process. Position-velocity diagrams were created 
by plotting the declination offset as a function of velocity chan- 
nel, where only coordinates with an uncertainty in the declina- 
tion of less than 0'.'5 were taken into account. Subsequently, the 
signatures of Keplerian rotation around a point source were fitted 
to the position-velocity diagrams (Fig.|7]i. In the fitting process, 
the values for M star (the mass of the central object), the central 
velocity of the system with respect to the local rest frame, and 
the declination offset from the source to the phase centre were 
varied, and a global minimum y 2 was determined. 

Assuming a Keplerian disk seen edge-on (i - 90°), a central 
mass of M star = 0.09 + 0.03 M G is found for IRS 63. The exact 
inclination of IRS 63 is hard to determine due to the high ex- 
tinction in the region. However, an inclination larger than ~ 45° 
can be ruled out because of the source's brightness at 3-5 /mi 
( Pontoppidan et al.|2003 i. Assuming an inclination of 30°, and 
taking into account the sin 2 i dependence, the mass of the cen- 
tral object increases to M star = 0.37 ± 0.13 M . For Elias 29, the 
points attributed to the emission from the dense ridge (Fig. |7]i 
were disregarded in the fitting procedure. A lower limit to the 
central mass of M stal = 0.62 + 0.14 M Q is found, under the as- 
sumption i = 90°. The flatness of the SED limits the inclination 
of this source to < 60° ( |Boogert et al.||2002"| ). A fully face-on 
orientation, on the other hand, is unlikely, given the presence of 
low surface brightness scattered K-band light (Zinnecker et al. 
1988 ). An inclination i = 30° yields a mass for the central source 
of 2.5 + 0.6 M Q . With such a high mass, the central star is likely 
to emerge as a Herbig Ae/Be star, once the surrounding envelope 
is dissipated. 

It is interesting to compare the inferred masses to those in- 
ferred for optically visible pre-main sequence stars. With a stel- 
lar mass of 0.37 Ms and a stellar luminosity of 0.79 L , IRS 63 
would be found to have an age of about 5 x 10 5 yr, when com- 
pared to evolutionary tracks of classical and weak-lined T Tauri 
stars in the p Ophiuchi cloud (Wilking et al. 2005 D'Antona 
|& Mazzitelli|| 1 997)> . A si milar age would be found for Elias 29 
(e.g., |Palla & Stahler|1993| l. 



5. Conclusions 

We used the SMA to observe the Class I YSOs Elias 29 and 
IRS 63 at 1.1 mm. The main results are as follows. 

- Both sources are detected in the continuum and in the HCO + 
J = 3-2 line, but not in the HCN J = 3-2 line. The HCO + 
emission toward IRS 63 is compact and associated with the 
YSO, whereas a significant part of the molecular emission 
in the direction of Elias 29 is due to an enhancement in the 
ridge from which Elias 29 formed. 

- Assuming that the continuum emission toward the sources is 
optically thin, and that any circumstellar envelope is resolved 
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Fig. 7. Velocity of the HCO + J = 3-2 line as a function of the 
declination offset from the phase centre for IRS 63 (upper panel) 
and Elias 29 (lower panel). The solid lines show the results from 
X 1 fits to the data, yielding central masses of 0.37 M Q (IRS 63) 
and 2.5 M Q (Elias 29) for inclinations of 30°; the dashed lines 
show formal 99% confidence intervals. Points attributed to the 
emission from the dense ridge are indicated and are not included 
in the fit. The declination offsets are different for the two panels. 



out by the interferometer on the longest baselines, we find 
empirical disk masses of 0.055 for IRS 63 and < 0.007 M Q 
for Elias 29. By comparing our interferometer data with 
single-dish observations, envelope masses of 0.058 M Q are 
derived for both sources, within 30" radius for IRS 63 and 
20" radius for Elias 29, where the value for Elias 29 should 
be treated as an upper limit because of the presence of the 
ridge enhancement. 

- Using a radiative-transfer code to create SEDs and ampli- 
tude vs. (u, v)-distance plots for different combinations of 
disk and envelope masses, we find disk masses of 0.13 M Q 
for IRS 63 and 0.004 M Q for Elias 29, and envelope masses 
of 0.022 M and 0.025 M . 

- The modelling results yield rather different ratios of 
Meny/Mdisk: 0.2 for IRS 63 and 6 for Elias 29, with uncertain- 
ties of a factor of a few. From this it is concluded that IRS 63 
is well on its way to become an optically visible Class II 
source, rather than an embedded Class I source. Elias 29, on 
the other hand, is more like the deeply embedded Class 
sources. 



- Velocity gradients in the HCO + J = 3-2 emission are in- 
terpreted as signs of Keplerian rotation, and indicate central 
masses of 2.5 ± 0.6 M for Elias 29 and 0.37 ± 0.13 M Q for 
IRS 63, for inclinations of 30°. These masses correspond to 
ages of a few 10 5 yr. 

With a larger set of sources treated in a similar manner as 
is presented for IRS 63 and Elias 29 here, it will be possible 
to constrain models for the evolution of the envelope, disk, and 
stellar masses through this important stage of YSOs. 
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